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Abstract 

We examine constraints on a simple neutrino model in which there 
are three massless and three massive Dirac neutrinos and in which the 
left handed neutrinos are linear combinations of doublet and singlet 
neutrinos. We examine constraints from direct decays into heavy neu- 
trinos, indirect effects on electroweak parameters, and flavor changing 
processes. We combine these constraints to examine the allowed mass 
range for the heavy neutrinos of each of the three generations. 
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1 Introduction 



Many models of neutrinos have been proposed to accomodate light or mass- 
less neutrinos. In a model with no right handed neutrinos, it is clear that 
neutrinos are massless. However, if there exist additional states which can 
play the role of Dirac partners to the left handed states, it is perplexing why 
neutrinos should be massless, or at least much lighter than their charged 
counterparts. Of course, neutrinos can be given small masses by coupling 
them to the standard Higgs doublet because of an extremely small Yukawa 
coupling, but it is more compelling to have an explanation for their small 
mass. A common explanation is the so called "see-saw" mechanism, in which 
the neutrinos remain light because the additional right handed states have a 
large Majorana mass. In such a model, neutrino masses are naturally small, 
since they are suppressed by the ratio of Dirac to Majorana masses, which 
is generally taken to be small. 

In this paper, we consider another viable alternative (see for ex. [2] , [3] , [4] ) . 
In addition to the three "right handed" neutrinos, there are three additional 
singlet particles. A lepton symmetry is imposed so that the only allowed mass 
terms are Dirac masses coupling the right handed neutrino to the standard 
left handed neutrino and to the additional singlet states. The conequence is 
that there are three heavy Dirac neutrinos, with mass determined primarily 
by the large mass term connecting the singlet and right handed neutrinos and 
three exactly massless neutrinos, the states orthogonal to the massive ones. 
Such a model has been considered before in several contexts; most recently 
it has been considered in the context of an Extended Technicolor Model with 
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a GIM mechanism [4]. In this type of model, the additional neutrino states 
could be quite light, on the order of a GeV. 

However, there are many constraints on such neutrinos. They are con- 
strained from direct searches for particles which have them in their final 
state, by universality constraints, and by flavor changing constraints. Cos- 
mological arguments are often used to constrain neutrino masses, but the 
neutrinos of this model are unstable so they do not apply. In this paper, we 
put these constraints together, making reasonable assumptions on the form 
of the mass matrix and mixing angles, to determine the allowed parameter 
regime. Many of these constraints apply quite generally to any model in 
which the left handed neutrinos mix with singlet states. Similar bounds were 
considered in ref. [18]. This paper updates the bounds, integrates them with 
those from LEP, and incorporates flavor changing bounds. We find with rea- 
sonable assumptions described below, the lightest neutrino can be as light 
as 2 GeV, although the third generation neutrino should be much heavier, 
greater than 80 GeV. 

The organization is as follows. We first present the model and we describe 
the approximations which we use to reduce the parameter space. We then 
consider constraints from meson and Z decays. Following this, we discuss the 
constraints from the fact that Gp will not have the same relation to standard 
model parameters when the muon cannot decay to the heavy neutrino state. 
We then look at flavor changing processes, which are in general permitted 
when no flavor symmetries are assumed. However, we assume mixing angles 
similar to those of the standard KM matrix, so there are approximate U(l) 
symmetries present. We then put together the constraints and consider three 
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models which describe the ratio of masses of the heavy neutrinos to determine 
the allowed parameter regime. Finally, we conclude. 

2 The Model and Simplifying Assumptions 

Many models have incorporated the neutrino scenario we discuss here. For 
example, it has been incorporated into GUT models [2], [3]. More recently, it 
has been shown how to incorporate such a model in an Extended Technicolor 
scenario [4]. We only consider the pheomonology of the lepton sector here, 
so we neglect the origin of the model and focus on the neutrinos. 

The standard model is extended by introducing three new left-handed 
neutrinos Sl and three right-handed neutrinos vr. Both left handed neu- 
trinos are coupled to the right handed neutrinos through Dirac matrices. 
All other possible mass entries are forbidden by a lepton number symmetry. 
Thus, 

-C mass = (u R 0)^ o)(^) +/? - c - 

This coupling results in three massive Dirac neutrinos and three massless 
eigenstates. The mass matrices D and S have different mass scales. The 
scale for D is constrained by SU(2) symmetry breaking whereas the scale for 
S is not, so it is reasonable to expect the masses in S to be larger. 

The mass of the heavy neutrinos is essentially determined by S. The 
electron, muon and tau neutrinos are a superposition of massless and massive 
eigenstates. The mixing to the massive neutrinos will, however, be small; 
it will be of the order of M D jM Sl where M D and M s are typical masses 
in D a To see more precisely how this mixing occurs, we need to find the 
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three massless eigenstates v° as well as the three with mass v . The mass 
matrix can be diagonalised by multiplying on the left and the right by unitary 
matrices: 

C X °y) * (£)-U(£). 

The unitary matrix V diagonalises DD^ + SS^ to give M 2 . The unitary 
matrix U is given by 

TT_ (UlA s V« ulA D V* \ 
^-KuIAdV* -UlAsV*)' 

where the matrices V, Up and Us are unitary matrices. They diagonalize 

M~ 1 VD and M~ 1 VS (where M~ l is the inverse of the diagonal mass matrix 

M) to give the diagonal matrices A D and A s : 

V^{M- 1 VD)U] } = A D 

V'^M^VSjUl = A s . 

The fact that the same V appears on the left for both these diagonaliza- 
tions is a consequence of the fact that the two matrix products M~ X V DDW*- M~ x 
and M~ l V SSW^ ' M~ x commute with each other, which follows in turn from 
the fact that their sum is the unit matrix. The most important part of U is 
the top right 3 by 3 block W D AsV'^ which links the electron, muon and tau 
neutrinos to the massive neutrinos u H . 

To extract bounds on the mass scales of S and D we need to make some 
simplifications to reduce the number of parameters. We will make the sim- 
plification that the matrices D and S are diagonalized by the same unitary 
matrices. In this case V = h X 3- If we then redefine the fields Sl by a 
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unitary transformation, absorbing the unitary matrix Us, we can rewrite the 
matrix U as: 



with A 2 D + A s = 7 3x3 . 

In this model the mass scale of the Dirac mass S is assumed to be much 
higher than the scale of the Dirac mass D. If this difference is sufficiently 
large, we can make the further simplification that As = 13x3- From here on 
the subscript Don A c will be dropped. 

At this point, we still have a large number of parameters. We simplify by 
assuming the matrix U D is similar in structure to the KM matrix for quarks. 
We notice that if there were no singlet left-handed neutrinos S L the matrix 
Ud would be the lepton equivalent of the KM matrix in the quark sector. 
We take the individual elements to be of the same magnitude as those in the 
KM matrix for quarks. 

We will use these approximations from now on. They leave six free param- 
eters: Mi, the masses of the heavy neutrinos and M D ., the masses induced 
by the mass matrix D which are defined as Mn i = Aj x Mj. In the following 
sections we will use experimental results to put limits on these masses. 

3 Direct Searches for Heavy Neutrinos 

Many searches for massive neutrinos have already been conducted. Massive 
neutrinos have been sought in the decays of n + [5] [6] [7], K + [8] [9] and 
charmed mesons [10] [11] [12], as well as in the neutral current production 
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of neutrino anti-neutrino pairs from e + e [13] [14] [15] collisions, and, more 
recently the decay of the Z [16] [17]. 

3.1 Meson Decays 

If it is kinematically allowed, any process involving the production of neutri- 
nos will be a source of heavy neutrinos. The creation process, however, will 
be supressed since the weak eigenstate neutrinos contain only a small mixing 
of the heavy neutrinos. Leptonic decays of mesons are thus one place to look 
for heavy neutrinos. 

At the lower end of the mass scale heavy neutrino creation in the decay 
of 7r + mesons has been investigated in references [5] [6] [7], and those of K + 
mesons in references [8] [9]. These experiments attempted to measure the 
mass of any heavy neutrino as it was created. This was achieved by stopping 
the 7T + and K + mesons and observing the energy of positrons emitted in 
their decay. The method did not rely on any assumptions about how the 
heavy neutrinos decayed. For massive neutrinos with masses less than 300 
Mev these experiments placed strict limits on the mixings, \U ei \ 2 and \U^\ 2 
of a heavy neutrino vf into the electron and muon neutrinos. For a range of 
masses the matrix elements |£/ e j| 2 and \U^i\ 2 were constrained to be less than 
10~ 5 . Since we assume the matrix Ud is almost diagonal, we can get direct 
bounds on Ai and A 2 of: Ai < 3 x 10~ 3 for the mass range 35 to 360 MeV, 
and A 2 < 3 x 10~ 3 for the mass range 80 to 325 MeV (see Figs(l) and (2)). 

Further limits on \U ei \ 2 and |t^| 2 come from the decays of charmed D 
mesons, [10] [11] [12], (see Figs(l) and (2)). Although similar searches can in 
principle be performed with decays of the B, they have not yet been done. 
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In total the limits from meson decays, with a few gaps, restrict Ai < 3 x 10~ 3 
for the mass range 35 MeV to 2 GeV, and A 2 < 3 x 1CT 3 for the mass range 
80 MeV to 2 GeV. For much of these ranges the bounds are much stricter 
than this. 

3.2 e + e Collisions at Low CM Energy 

Heavy neutrino anti-neutrino pairs would be created by weak interaction 
currents in e + e~ annihilations, but since the the center of mass energy of 
these collisions is less than the W and Z mass the cross section for the 
creation of these neutrinos is extremely small. Experiments [13] [14] [15] were 
aimed at detecting the decay of a heavy fourth generation neutrino and they 
thus made the assumption that the heavy neutrino had the same coupling 
to the Z and W as the other neutrinos. This is not the case for the model 
studied in this paper where each heavy neutrino introduces a mixing angle 
factor of \Uu\ 2 into the weak interaction couplings. Reinterpreting the data 
of these experiments including the extra mixing angles results in constraints 
that are negligible in comparison to the other bounds studied in this paper. 

3.3 Z decays 

Massive neutrinos, lighter than M z , would also be created in Z decays, and 
the experiments [16], [17] have already conducted searches for heavy isosinglet 
neutrinos; the type discussed in this paper. The most abundant supply of 
heavy neutrinos would come from the decay of a Z into one heavy neutrino 
and one massless neutrino. For a Z decaying into a heavy neutrino vf (lighter 
than the Z) and any of the three massless anti-neutrinos z/° the creation is 
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suppressed by: 



j'=l l=e,fi,r 



where R n has the following meaning: if N is the number of neutrinos (from 
one family of the standard model) created in the experiment then then num- 
ber of heavy neutrinos, vf , created is RiN. 

Experiments aim to detect the neutrino by its decay. The decay of the 
neutrino would be quite distinctive. In general, it will decay to a high energy 
lepton and a virtual W or Z, which would then decay into leptons, or hadrons 
if the neutrino is massive enough. The total decay rate can be written in 
terms of the rate for muon decay as follows [18]: 



T(z/f -> leptons /hadrons) = E U ^+3)| 2 (tt) ®i( M i) r (v ~> evv ) 

where $;(Mj) is a factor that weights the decay rate for a single channel by 
the effective number of channels into which there is sufficient energy to decay 
and takes into account the different Feynman diagrams. 

There are two reasons why decays like this might not have been seen in 
experiments: 

Very few heavy neutrinos are produced. If we assume that nearly all 
the neutrinos decay inside the detector, then the fraction, Rj, that 
decay inside the detector is given by: 

Rd = R% = A 2 . 

If Ri is sufficiently small, no neutrinos would be detected. 
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The neutrinos have a long lifetime. If the neutrinos are light, they could 
have a very long lifetime and thus decay almost entirely outside of the 
detector. 

We can then calculate the fraction, R d of Z's that would decay inside the 
detector: 



where 7 is the time dilation factor due to the relativistic motion of the neu- 
trino, and for M z » Mi is given by 7 = M z /2il /,. Sd is the size of the 
detector and c is the speed of light. 

The experiment of reference [16] involved a search through 4x 10 5 hadronic 
Z decays and placed limits of Aj < 0.014 over the range 5 to 50 GeV. Above 
50 GeV the phase space for heavy neutrino production becomes smaller and 
the limits placed on the Aj become less strict. Below 5 GeV the limits are 
reduced due to the long lifetime of the neutrinos. In Figs (1), (2) and (3), 
are marked out the forbidden regions in the M i: ^- plane for the three heavy 
neutrinos. 

4 Changes In Weak Interaction Decays and 
Parameters 

Aside from direct searches for the heavy neutrino, the existence of the heavy 
neutrino will affect precision measurements of various electroweak processes. 
This can be the case because G F will no longer have the standard model 
relation to sin6w, since the muon decay rate will be different if the muon 
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cannot decay into the heavy neutrino. This would change the relation be- 
tween precisely measured electroweak parameters, for example the W mass 
or sin 2 9 W as measured in the forward-backward asymmetry. Furthermore, 
it would lead to an apparently nonunitary KM matrix. 

Further constraints come from pion decay branching ratios if the heavy 
electron or muon neutrinos are heavier than the pion. Similarly, universality 
could be violated and would be seen in tau decay. Finally, the Z width can be 
affected, both indirectly though a change in the extracted sin 2 9, and directly 
if the neutrinos are heavier than the Z. 

4.1 Muon decays and the Fermi coupling constant Gf 

The Fermi coupling constant Gp is the effective coupling constant for four 
fermi interactions and is measured extremely accurately from muon decays. 
If the mass of the massive neutrinos is greater than that of the muon, the 
decay width for the muon would be decreased, since it would not be able to 
decay into the heavy neutrinos; this in turn would lead to a change in the 
predicted value of Gp. 
Specifically: 

3 

which leads to 

K G FTmuon decays — + ^1) 

where 5 means the fractional change. Of course Gp is a measured num- 
ber. What we mean here is the change in the coefficient of the four fermion 
operator which yields muon decay. 
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4.2 Semileptonic Decays and the KM Matrix 

The estimates of the semileptonic processes would also be affected but to 
a lesser extent. The same value of Gp is also used for the effective cou- 
pling constant for semileptonic decays, where elements of the KM matrix are 
determined. One would expect these elements to be part of a unitary matrix. 

The important point to consider is that the effective coupling constants 
for the leptonic and semileptonic four fermion interactions would no longer be 
the same, and if it was assumed that they were, the predicted matrix elements 
for the KM matrix would no longer be those of a unitary matrix. We can 
check the unitarity of the KM matrix by looking at the matrix elements 
(KM) ud , (KM) US and (KM) ub ; the sum of their square magnitudes must 
add up to one. The effect of having heavy neutrinos would be to make 
this sum slightly bigger than one. The most important shift will come from 
the change in nuclear beta decays used to determine the (KM) ud element. 
Consequently, what must be compared are the changes in the value of Gp 
and in the rates for nuclear beta decays. 

Specifically, as above, 

3 

(GF) new = \UeiU jfl \ (G F ) old , 

which leads to 

KG Ffrauon decays ^ ~(A? + k\) 

and similarly 

5 (beta decay) ~ — 

where 5 means the fractional change. The fractional change in the width of 
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the muon minus the fractional change in nuclear beta decays must be less 
than the experimental uncertainty in the sum of the matrix elements. From 
reference [19]: 

\(KM) ud \ 2 + \(KM) US \ 2 + \(KM) ub \ 2 = 1 (+8.6 x 10" 4 , -4.7 x 1(T 3 ). 
This leads to A 2 < 6 x 1CT 2 (2a) if the neutrinos are heavier than the muon. 
4.3 Mw and sinOw 

Changes in Gp would also affect the prediction of other weak interaction 
parameters. The ratio of the mass of the W and Z for example depends 
upon Gf- Specifically [19]: 



where 5v is a radiative correction parameter much less than one. Using the 
above the change in the predicted value of Mw/Mz due to the change in Gf 
is: 



Current experimental bounds [20] place 5(Mw/Mz) < 7.7 x 10 3 (2c) which 
gives a bound for Gf of: 



In fact this bound is too strong due to the uncertainty in the top quark mass. 
However, since it iss less strict than the bound coming from the KM matrix, 
it will not be incorporated. 




8(M W /M Z ) = 0.088 x S(G F ). 



S(G F ) < 8.8 x 10" 2 (2a) 
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The Weinberg angle sin 2 9w also depends on (the on shell definition is: 
sin 2 9 W = (1—Myy/Mz) and this can be compared with the forward backward 
assymetry of the process e + e~ — > //, which depends on 6 W . However, this 
too is weaker than the constraint from the unitarity of the KM matrix. 

4.4 Pion Decay Branching Ratios 

The ratio of the two decay channels for a n ± ; n — > e v e and ir — > /i provides 
another bound [21]. In this model: 

ll n " eVe \ = 1-233 x 10- 4 (1 - f(l^| 2 " l^l 2 )), 

where the factor 1.233 x 10~ 4 is the theoretical value of the ratio in the 
standard model[22]. The mixing angle factors apply for neutrinos too heavy 
to be produced. Experimentally the ratio is known to be: (1.218 ± 0.014) x 
10~ 4 [23]. Using the fact that U D is almost diagonal and that A 2 < 6 x 10~ 2 
leads to: 

0.035 > (A 2 - A 2 ) =^ Ai < 0.18 (2a) 

4.5 Tau decays 

If the neutrinos are all heavier than the tau then the decay width of the 
tau would also be affected. As for the case of the muon decay it can be 
shown that the partial width T(r — > ev T v e ) would be reduced by ~ (A 2 + A§) 
and the partial width T(r — > iiv T v mv ) would be reduced by ~ (A 2 , + A3). 
Consequently, the partial width for decay into leptons would be reduced by: 

~4a 2 + 1a 2 + a 2 ). 

13 



This fractional change minus the fractional change in the width for muon 
decay must be less than the experimental uncertainty of the partial width 
for the tau. This gives a further bound on the A,: 

||A? + |A1- A§|< 0.015 (la) 

where the uncertainty in the partial width of the tau is 1.5% [19]. To obtain 
a bound for A3 we use the following formula for calculating the error at the 
la level of a sum of terms each with their own errors: 

A^ = ((^)L + (^)L + 0.015f , 

where, at the la level, we use the bounds from the previous section Af < 
0.017 and A 2 2 < 0.0018. This leads to a bound on A 3 at the 2a level of: 

A 3 < 0.18 (2a), 
4.6 The Width of the Z 

For neutrinos heavier than the Z the width will be reduced, since the decay 
into the heavy neutrinos will no longer be kinematically allowed. Experimen- 
tally, the partial width, Y v9 of the Z is known to an accuracy of 1.8% [24]. In 
this model 

r,poci£|i-A, 2 l = i-iE^, 

i i 

where the sum over % is only over neutrinos heavier than the Z. This gives 
the bound: 

]TA, 2 < 0.108 (2a). 
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There is also an effect if muon decay is changed. However it is not numerically 
as important. 

In Figs (1), (2) and (3) are plots of the bounds placed on the Aj by all 
the processes considered in sections four and five. 

5 Lepton Flavor Changing Processes 

Flavor changing processes were also examined in this model. These processes 
are exactly analagous to flavor changing processes in the quark sector. Three 
processes with strong experimental bounds were considered [19]: 

(i) /i -> e 7 ; experimentally: ^^gy < 5 x KT 11 , 

(ii) n — > ee + e~; experimentally: r ^^^!- ) - < 1CT 12 , and 

(iii) pTi ^eTi- experimentally: r ^^ e) < 5 * 10" 12 . 

These processes can only occur via loop diagrams involving the exchange of 
virtual neutrinos. The couplings of the neutrinos to the muon and the elec- 
tron involve the unitary matrix U; specifically, the neutrino- W- muon vertex 
includes a factor Uj M for coupling to the ith neutrino, and the factor \] ei is 
included with the neutrino-W^-electron vertex. The amplitudes are obtained 
by summing over all intermediate states i. All terms proportional to the 
sum Z)i=iU e jUl M are automatically cancelled since U is unitary. This is the 
GIM mechanism. It is an analog of the strong suppression of neutral current 
flavor changing processes in the quark sector. Notice this is independent of 
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the approximations we made. 

In all cases there is very strong GIM suppression. Flavor changing pro- 
cesses can only proceed via an intermediate heavy neutrino state. But the 
coupling to the charged neutrinos is then suppressed. Therefore these con- 
straints will only dominate in the region where there are no other strong 
constraints, namely for neutrinos more massive than the Z . 

For the purpose of calculations, the masses of the electron, muon and tau 
and the mass matrix D are generated in the standard way by coupling to 
the Higgs, so that the loops involved charged Higgs. t'Hooft gauge is used 
throughout, simplifying the form of the propagators and setting the masses 
of the W and the charged Higgs to be the same. We will now consider, in 
detail, the three flavor changing processes. 

5.1 fi — > e7 

A calculation for the case where the neutrino masses are much less than Mw 
is elaborated in Cheng and Lee [25]. Below we will give an outline of the 
calculation for the general case, where the neutrino masses are not assumed 
to be less than M w . To simplify the calculations the electron is taken to 
be massless. The general form of the amplitude is constrained by gauge 
invariance, thus the gauge invariant form of the amplitude with a massless 
electron is given by 

< e,i\(S- >= Au e [(l - ^)ik v e x a Xu i0\u^ 

where k is the photon four momentum, e is the polarization of the photon and 
A is a constant to be determined. The partial derivative term is included to 
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ensure that, in the hypothetical case where the muon mass goes to zero, only 
the lefthanded component of the muon coupled to the W is involved in the 
decay. The amplitude can then be rewritten using the Gordon decompostion 
as: 

< e, 7| (S - 1) \n >= A m^u e [(l + 75) (2e.p - m^\u^ 

where m M is the muon mass and p is the four momentum of the incoming 
muon. To simplify the calculation, only the terms proportional to e.p need 
to be calculated. In principle there are 8 possible diagrams contributing (see 
Fig. (8a)). Diagrams 4 to 8 contain only the j. term and thus can be ignored. 
They will cancel with similar terms coming from the first four diagrams. 

In evaluating diagrams 1 to 4 we can define a factor i] for each of the 
diagrams i =1 to 4, and for each of the neutrinos j —1 to 6. Summing over 
the six neutrinos j the contribution of diagram i to the constant A is 

-iK^U^Uliy^M^) 

3=1 

where K = — e 3 /(2 sin 2 9w)- We can also define the sum 

1=1 j=i 

so that the total contribution to the constant A from all the diagrams is 

-^ e /(32vr 2 M^). 

Performing the calculations gives the following results for the i] : 

I] = -[ aj {a} - 3a, + g) + a^aj - |)o 2 logo,], 
Ij = ~[\a) - |a 2 - j^aj + I + \a){aj + l)S]logSj], 
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If = 0, and 

where the variable 8j is given by Sj — for j — 1, 2, 3 and by Sj = (M U( j_ 3 ^/M w ) 2 
for j = 4, 5, 6, and aj is defined as aj = 1/(1 — 5j). Using the expression for 
U, the sum 7 Me is then given by: 

6 

I fie = H U D e(j - 3) U D( ._ 3)fi A U - 3) 
j=4 



X [ " 2 3 + T°? " T % + 2 " 2°^ ° g ^ 



This can be approximated for the two cases where the neutrino masses are 
all either much less than or much greater than the mass of the W. Thus, 

he * - Ej=i U Dej ul.An 4 S {j+3) if all Mj < M w 

* £j=i^<A?(5 + ^T + §kg<W if all M,- >M W . 

Averaging over the initial spins and summing over the final spins and mo- 
menta leads to the decay rate: 

_ aG 2 F ml 
^ n 128tt 4 1 Vl ' 

where ct is the fine structure constant. This can be compared to the decay rate 

G 2 m 5 

T^e7 = 1 g 27r ^ to get the ratio which can then be compared to experiment 
to get an upper bound on I^ e : 

T(/i -> e 7 ) 3a 2 in _i 



< 1(T 1U (2a), 



r(/x — ► ez/z/) 27r 

=}► |/ Me | < 1-7 x 1(T 4 (2a). 
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5.2 fi — > ee + e 

This process can occur via extensions of the diagrams in \i — > where the 
7 is virtual and splits into an electron positron pair, or it can take place 
via box diagrams (see Fig (8b)). The contribution from box diagram(l) 
dominates the other box diagrams. This is due to the fact that the other 
diagrams involve exchange of virtual charged Higgs whose coupling to the 
electron and muon is supressed by a factor of the order of the Dirac mass D 
over My/. The contribution from the extensions of the diagrams in fi — > &-/ 
is of the same order as the first box diagram, but for an order of magnitude 
estimate we approximated the whole amplitude from the first box diagram 
only, (diagram(l) Fig(8b)). 

The amplitude for the box diagram is calculated using the approximations 
that the neutrino masses are much greater than the muon mass, and that 
the electron mass is zero. Using these approximations the amplitude for the 
process is: 

• 4 

A = "32tt2 sin 4 Q wM 2 v ^ P ~W P u q2 P-l v v q3 J, e , 
where the dimensionless factor J^ e is: 

6 

J^e = ]T UejU^Al^aj - 1 + a)5j log 5,) 

3=1 

with the variable 5j given by 5j = for j = 1, 2, 3 and by 5j = (M U( j_ 3 ^/M w ) 2 
for j = 4, 5, 6, and aj defined as dj — 1/(1 — 5j). As for the previous flavor 
changing process, we can make approximations for the cases where all the 
neutrino masses are less than the mass of the W and where they are all much 
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greater than the mass of the W. We obtain: 

Jue ^ E?=i Uo^Ul.^S^il + log^+s)) if all Mj < M w 
* E?=i U De3 tf Dj A] - [1 + jj-Ml - log% +3) )] if allM, > Mwr. 



Averaging 2 over the initial spin states and summing over the final spins 
and momentums leads to the decay rate: 

Gpin^a 2 2 
r ^ ee+e " = 7687r5sin 4 ^ |JMel ' 

where a is the fine structure constant. As before this can be compared to the 
decay rate T^ euP = to obtain the ratio which can then be compared 

to experiment to get an upper bound on J^ e : 



r(/x — > ee + e ) a 2 
T(fj, — > evv) Att 2 sin 4 6 W 



J^Y < 2 x 10~ i2 (2a), 



| J Me | < 2.8 x KT 4 (2a) 



5.3 /iTi^eTi 



O.Shanker [26] has performed some careful calculations for fie conversion for 
different nuclei. These calculations involve using an effective Hamiltonian for 
the muon-electron-q-q vertex where the q's represent either two up quarks or 
two down quarks. This effective Hamiltonian is obtained from box diagrams 
very similar to the ones in the previous section except that the outgoing 
electrons, labeled with momenta q2 and q^, are replaced by an incoming and 
outgoing quark in the Titanium nucleus, see Fig(8c). The calculation for the 



20 



amplitude from the previous section can be carried over with very few changes 
to give the amplitude for the process involving the up quark; fxu — > eu: 

G F aJ^ e x - p - 

V 24:71 sin O w 

and the same amplitude for the down quark. These calculations assume that 
the quark mass is less than the neutrino mass. We can then write an effective 
Hamiltonian for this interaction: 

2G 

H e ff = g^^eP~'j X fi(uP~'j x u + dP'^d), 

where g = aJ (Ue /(47rsin 2 9 W ). Using the calculations of O.Shanker [26], we 
can then obtain the ratio between the decay rate for \x Ti — > e Ti to the rate 
for muon capture by the nucleus which can be compared to experiment to 
obtain another bound on J Me : 



2 



265.64 . . |J Me | 2 < 10" n (2a 



T(fi Ti capture) 167T 2 sin 4 6^ Me 

|J Me | < 7.6 x (2a). 

The diagrams for the two processes, uTi —> eTi and muon capture, are 
essentially the same as the diagrams for u — > ee + e~ and u — > euu but the 
ratio of the decay rates of the first two processes is much greater than for the 
second two. Thus although the bounds placed by experiment on uTi — > eTi 
are not as strong as those for u — > ee + e~ , it is the process uTi — > eTi which 
places the strongest bounds on the size of J Me . 

This difference can be explained by coherence effects. The dominant 
process for \xTi — > eTi leaves the Ti nucleus in it's ground state [26], which 
is a coherent process involving summing the amplitude over all the nucleons. 
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Muon capture, on the other hand, is an incoherent process involving summing 
the square of the amplitude over all the protons. 

6 Discussion and Conclusions 

One of the desired results of this model was that it would provide a scenario 
in which weak interaction symmetry breaking could give the neutrinos a mass 
matrix on a scale similar to that of the electron, muon and tau, while still 
maintaining massless neutrinos. To investigate this all the plots discussed 
in this section are marked with a dashed line along which the masses M D . 
induced by the mass matrix D (from weak interaction symmetry breaking) 
are the same as the electron, muon and tau. For all the plots the excluded 
regions lie to the left of the curves. 

The bounds from sections three and four are plotted separately for each 
of the Ai in Figs (1), (2) and (3). To fulfill the scenario in which M Dl = M e , 
M D2 = and M Ds = M T we see that the mass of the third neutrino must 
be greater than M W) the second must be heavier than 10 GeV and the first 
heavier than 2 GeV. 

To examine if there are further restrictions from the flavor changing pro- 
cesses of section five, the Aj have to be plotted on the same graph since the 
factors I^e and J^ e are functions of all three Aj. In fact due to the very small 
mixing of the third massive neutrino vf into the electron and tau neutrinos 
(the mixing matrix Ud is chosen in this analysis to be like the KM matrix) 
I^e and J Me are virtually independant of A 3 . 

Figs (4), (5), and (6) plot out the constraints from the flavor changing 
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processes for three different scenarios. They all assumed that the masses 
M D . generated by the mass matrix D were in the same ratio as the masses 
of the electron muon and tau ie. M Dl : M Da : M 3 = M e : M M : M T . The 
dotted line, as before, marks out the line along which the masses M D . are 
actually the same as the electron, muon and tau masses. The flavor changing 
processes are plotted alongside all the other constraints from sections three 
and four. It is immediately clear that flavor changing processes do not rule 
out any of the line along which M Dl = M e , M D2 = and M Ds = M T . The 
bounds from Z decays and, for the lightest neutrino, meson decays are much 
more important. 

6.1 Three scenarios 

In the scenarios that follow four different ratios of the neutrino masses M { are 
considered. The bounds given at the end of the discussion of each scenario 
assume that M Dl = M e , M Di = M M and M 3 = M T . 

Scenario 1 Fig(4) .U, M 2 M 3 . 

In Fig(4) are plots of the allowed regions taking into account all the 
experimental constraints from sections three, four and five. Areas to 
the left of the curves are ruled out. The plot is of the mass M 2 of the 
second heavy neutrino against ^, the ratio between the two mass scales 
generated by S and D. The most important constraint comes from the 
limits set by Z decays on the third neutrino. If they are to lie on the 
dashed line all three neutrino masses are constrained to be greater than 
M w . 



23 



Scenario 2 Fig(5) M 1 : M 2 : M 3 = 1 : 15 : 60. 

Again, areas to the left are ruled out by experiment and the plot is of 
the mass M 2 of the heaviest neutrino against the ratio between the 
masses M 2 and M D2 . In this case the most important constraints are 
those set by D decays on the mass of the first neutrino and those set by 
Z decays on the mass of the third neutrino. If the masses are to lie on 
the dashed line M 2 must be greater than 30 GeV. Dividing this by 15 
and multiplying by 4 gives the bounds for the first and third neutrinos 
respectively. The bounds for the three neutrinos are thus: M 1 > 2 GeV 
which is equivalent to M 2 > 30GeV and M 3 > 120 GeV. 

Scenario 3 Fig(6) M : M> : M 3 M : M,, : M T 

In this scenario it is the constraints set by D decays on the mass of the 
first neutrino that are most important and the corresponding bounds 
for the three masses are (for masses lying on the dashed line): Mi > 
2 GeV, M 2 > 400 GeV and M 3 > 3500Ge\/. 

6.2 Conclusions 

In this paper we have examined the experimental consequences of a model of 
massive neutrinos and have excuded a large region of the parameter space. 
Specifically we have found that, in the scenario where the mass contributions, 
Mr>i, from weak interaction symmetry breaking are the same as those for the 
electron, muon and tau, the neutrino masses are approximately constrained 
as follows: 

Mi > 2 GeV, M 2 > 10 GeV, and M 3 > 80 GeV. 



24 



This means that either the Dirac mass connecting standard left handed neu- 
trinos to right handed neutrinos has entries less than their charged counter- 
parts, or one would not expect all the neutrinos to be light. It is clearly 
nonetheless of interest to improve the bounds. Clearly improved statistics at 
LEP will give stronger constraints. Furthermore, the bound on Mi can be 
improved by looking for heavy neutrinos in B decays. 
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Figures 



Figures 1, 2 and 3 show bounds placed on the heavy neutrino masses by 
considering all the constraints in chapters three and four. The plots are 
of the neutrino masses M Di against the ratio l/A Di = Mi/M Di . Any 
region to the left of a solid line is forbidden by experiment. The dashed 
line is the line along which the masses M D . are equal to the electron, 
muon and tau masses. 

Fig 1 Regions excluded(to the left of curves) in the Mi, plane from: 
(a)ref[5] Massive neutrinos in pion decays; (b)ref[9] and (c)ref[8] Mas- 
sive neutrinos in Kaon decays; (d)ref[10] Massive neutrinos in D meson 
decays; (e)ref[21] /section (4.2) Banching ratio r(7T — > eu)/Y(j[ — ■> jiu); 
(f)ref[16] section (3.4) Z decays. 

Fig 2 Regions excluded(to the left of curves) in the M 2 , ^ plane from: 
(a)ref[9] Massive neutrinos in Kaon decays; (b)ref[10] and (c)ref[10] 
Massive neutrinos in D meson decays; (d)Section(4.1.1) Changes in 
G F , (e)ref[16] section (3.4) Z decays. 

Fig 3 Regions excluded (to the left of curves) in the M 3 , ^ plane 
from: (a) section(4.1.2) Tau decays; (b)ref[16] section (3.4) Z decays. 

Figures 4, 5, 6 and 7 Regions excluded (to the left of curves) in the M 2 , 
plane from the experimental constraints from chapters three four 
and five. Each plot has the same ratio between the masses M Dl , M D2 
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and Md 3 induced by D, This is chosen to be Mo 1 : Mr> 2 : Md 3 = 
M e : M M : M T . The dashed line correponds to the line along which 
M Dl = M e , M D , 2 = M M and M Dz = M T . The different plots correspond 
to three different ratios of M x : M 2 : M 3 (the final neutrino masses). 

Fig 4 Case (1): M. M, M 3 . 

Regions excluded(to the left of curves) in the M 2 , ^ plane from: (a) All 
the restrictions in the Mi plane studied in chapters three and four; 
(b)All the restrictions in the M 2 ^ plane studied in chapters three 
and four; (c)All the restrictions in the M 3 plane studied in chapters 
three and four; (d)Bounds from /i — > e-y; (e)Bounds from /iTi — > eTi. 

Fig 5 Case (2): M x : M 2 = 1 : 15 : 60. 

Regions excluded(to the left of curves) in the M 2 , ^ plane from: (a) All 
the restrictions in the Mi plane studied in chapters three and four; 
(b)All the restrictions in the M 2 plane studied in chapters three 
and four; (c)All the restrictions in the M 3 plane studied in chapters 
three and four; (d)Bounds from \x — > e^y; (e)Bounds from \xTi — > eTi. 

Fig 6 Case (3): M ; : M- 2 : M :i M, : .\/ /( : M r . 

Regions excluded(to the left of curves) in the M 2 , ^ plane from: (a) All 
the restrictions in the Mi plane studied in chapters three and four; 
(b)All the restrictions in the M 2 plane studied in chapters three 
and four; (c)All the restrictions in the M 3 plane studied in chapters 
three and four; (d)Bounds from \i — > cy; (e)Bounds from fj,Ti — > eTi. 
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7 Feynman diagrams for the flavour changing processes of chapter five. 
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